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We discuss the final results of the MACRO experiment on atmospheric neutrino oscillations. The data concern 
event topologies with average neutrino energies of ~ 3 and ~ 50 GeV. Multiple Coulomb Scattering of the high 
energy muons was used to estimate the neutrino energy event by event. The angular distributions, the L/E^ 
distribution, the particle ratios and the absolute fluxes all favor — > v T oscillations with maximal mixing and 
Am 2 ~ 0.0023 eV 2 . Emphasis is given to measured ratios which are not affected by Monte Carlo (MC) absolute 
normalization; a discussion is made on MC uncertainties. A preliminary search for possible Lorentz invariance 
violation contributions to atmospheric neutrino oscillations is presented and discussed. 



1. Introduction 

MACRO was a large area multipurpose under- 
ground detector pQ designed to search for rare 
events and rare phenomena in the penetrating 
cosmic radiation [2]. It was located in Hall B of 
the Gran Sasso Lab at an average rock overbur- 
den of 3700 m.w.e.; it started data taking with 
part of the apparatus in 1989; it was completed 
in 1995 and run in its final configuration until 
the end of 2000. The detector had dimensions of 
76.6 x 12 x 9.3 m 3 and provided a total acceptance 
to an isotropic flux of particles of ~ 10, 000 m 2 sr; 
vertically it was divided into a lower part, which 
contained 10 horizontal layers of streamer tubes, 
7 of rock absorbers and 2 layers of liquid scin- 
tillators, and an upper part which contained the 
electronics and was covered by 1 layer of scintil- 
lators and 4 layers of streamer tubes. The sides 
were covered with 1 vertical layer of scintillators 
and 6 of limited streamer tubes. 

MACRO detected upgoing v^s via charged 
current interactions, — > //; upgoing muons 
were identified with the streamer tube system (for 
tracking) and the scintillator system (for time-of- 
fliglit) . The events measured and expected for the 
3 measured topologies, deviate from Monte Carlo 



*see Ref. [1] for a fist of MACRO Authors and Institutions 



(MC) expectations without oscillations; these de- 
viations and the L/E v distribution point to the 
same v T oscillation scenario [21 — [31 - Here we 

also present the results of a preliminary search 
for possible Lorentz invariance violation contri- 
butions to atmospheric neutrino oscillations. 



2. Atmospheric neutrinos. Monte Carlos 

The measured upthroughgoing muon data of 
Fig. were compared with different MC sim- 
ulations. In the past we used the v flux com- 
puted by the Bartol96 group The system- 
atic uncertainty in the predicted flux was esti- 
mated at ±17 %; this is mainly a scale error that 
does not change the shape of the angular distri- 
bution. A similar MC (Honda95) 9] was used by 
the SupcrK Collaboration (SK) (Hfl. Recently 
new improved MC predictions for neutrino fluxes 
were made available by the Honda (HKKM01) 
[TT] and FLUKA [T2J groups. They include three 
dimensional calculations of hadron production 
and decay and of neutrino interactions, improved 
hadronic model and new fits of the primary cos- 
mic ray flux. The two MC yield predictions for 
the non oscillated and oscillated fluxes equal 
to within few % 0. The shapes of the angu- 
lar distributions for oscillated and non oscillated 
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Bartol96, new FLUKA and new Honda fluxes are 
the same to within few %. The absolute values 
of our data are higher than those predicted by 
the new FLUKA and Honda MC, Fig. [TJ A sim- 
ilar situation is found in the new SK data |13j . 
The L3C and BESS cosmic ray results presented 
at this workshop [H] lead to v fluxes in better 
agreement with Bartol96 and Honda95 predic- 
tions. The evidence for neutrino oscillations is 
mainly due to the shape of the angular distri- 
bution and this is the same in all MCs. Also 
the ratios of the medium to high energy measu- 
rements and of the two low energy data samples 
are MC independent. Our data suggest that the 
FLUKA normalization should be raised by 25% 
at E v ~ 50 GeV and by 12% at E v ~ 3 GeV. 

3. MACRO results on atmospheric i/'s 

The upthroughgoing muons come from in- 
teractions in the rock below the detector; muons 
with > 1 GeV cross the whole detector. The 
corresponding v^s have a median energy of 50 
GeV. Many systematic effects and backgrounds 
were studied 0ZI The data, Fig. [TJi, agree m 
shape and absolute value with the oscillated Bar- 
tol96 MC, for Am 2 = 0.0023 eV 2 . 

— » v T versus — ► v s . The ratio 
Ri = Vertical/Horizontal = N(—l < cosO < 
-0.7) /iV (-0.4 < cos9 < 0) was used to test the 
— * v s oscillation hypothesis versus — > v r 2 
[S] [7|. The — > v s oscillations (with any mix- 
ing) are excluded at 99.8% c.l. with respect to 
— > v r oscillations with maximal mixing jjj. 
Oscillation probability as a function of the ratio 
L/E u . E v was estimated by measuring the muon 
energy, E^, by means of muon Multiple Coulomb 
Scattering (MCS) in the rock absorbers in the 
lower MACRO. The space resolution achieved 
was ~ 3 mm. The distribution of the ratio R = 
{Data/ MCnoosc) obtained by this analysis is plot- 
ted in Fig. El versus (L/E v ) 0. 

The Internal Upgoing (IU) muons come from 
~ 3 GeV z/^'s interacting in the lower apparatus. 
Compared to the no-oscillation prediction there 
is a reduction in the flux, without distortion in 
the zenith distribution shape, Fig. QJ). The MC 
predictions for no oscillations in Fig. ^ is the 
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Figure 1. Zenith distributions for MACRO data 
(black points) for (a) upthroughgoing (top) and (b) 
semicontained muons. In (a) there are comparisons 
between the predictions of old and new MCs with 
and without oscillations. In (b) the dashed line is the 
no-oscillation Bartol96 MC (with an error band) ; the 
solid line is for — » v T oscillations with sin 2 2i? = 1 
and Am 2 = 2.3 • KT 3 eV 2 . 



dashed line with a 21 % systematic band. 

The Upstopping (UGS) muons, due to ~ 3 GeV 
f M 's interacting below the detector, yield upgoing 
muons stopping in the detector. The Semicon- 
tained Downgoing (ID) fi's are due to f^-induced 
downgoing fx's with vertex in the lower MACRO. 
The 2 types of events (not shown in Fig. ^) are 
identified by topological criteria. 

4. Oscillation parameters. 

In the past, in order to determine the oscil- 
lation parameters, we made fits to the shape of 
the upthroughgoing muon zenith distribution and 
to the absolute flux compared to the Bartol96 
MC. The other data were only used to verify 
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Figure 2. Ratio (Data/MC Bartol96) versus the es- 
timated L/E v for the upthroughgoing muons (black 
circles) and the semicontained up-p (open circle). 
The horizontal dashed line at Data/MC=l is the ex- 
pectation for no oscillations. 



1Cr 1 - 



. 10- 



1 cr- 



ier' 



0.2 



0.4 



0.6 



0.8 



1.0 



Figure 3. Interpolated qualitative 90% C.L. contour 
plots of the allowed regions for the MACRO data us- 
ing only the ratios Ri,R.2, i?3 (outer continuous line) 
and using also the absolute values assuming the va- 
lidity of the Bartol96 fluxes (dotted line). 



the consistency and make checks. The result was 
Am 2 = 0.0025 eV 2 and maximal mixing j^j [3]. 
Later, also the L/E v distribution was used UJ. 

In order to reduce the effects of systematic un- 
certainties in the MC we recently used the follow- 
ing three independent ratios \7\ and we checked 
that FLUKA, Honda and Bartol96 MC simula- 
tions yield the same predictions to within ~ 5%. 

(i) High Energy Data: zenith distribution ra- 
tio: Ri = N vert /N hor 

(ii) High Energy Data, neutrino energy measu- 
rement ratio: R 2 = Ni ow /N hlgh 

(iii) Low Energy Data: 

ratio R 3 = {Data/ MC) W /{Data/ MC) ID+UGS . 

With these ratios, the no oscillation hypothesis 
has a probability P ~ 3 • 10~ 7 and is ruled out 
by ~ 5(7. By fitting the 3 ratios to the — > 
v T oscillation formulae we obtain sin 2 2i? = 1, 
Am 2 = 2.3 • 10" 3 eV 2 and the allowed region 
indicated by the solid line in Fig. [3J 

If we use the Bartol96 flux we may add the 
information on the absolute fluxes of the 

(iv) high energy data (systematic scale error of 

£17%) i?4 = N meas /NMCBartol- 

(v) low energy semicontained muons, with a 



systematic scale error of 21%, R 5 ~ 

Nmeas / N MCBartol- 

These informations reduce the area of the allowed 
region (dashed line in Fig. |3J), do not change the 
best fit values and bring the significance to ~ 6<r. 

We recall that in the 1984 proposal it was 
stressed that the new region in (Am 2 , sin 2 2(9) 
that MACRO would cover is as indicated in Fig. 

where the shape of the angular distribution of 
upthroughgoing /x's is very sensitive. 
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Figure 4. From the 1984 MACRO proposal: 1984 
limits on Am 2 , sin 2 2$. The shaded region indicated 
the expected improvement in the explored region ob- 
tainable with the MACRO experiment. 
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Figure 5. 90%C.L. upper limits on the LIV param- 
eters Av and sin 2 2$„, assuming Am 2 = 0.0023 eV 2 , 
sin 2 2$ m =1 for the mass induced oscillations (solid 
line). The upper (lower) dotted line refers to Am 2 = 
0.0015 (0.0034) eV 2 . The dashed line shows the limit 
obtained using the selection criteria of ref. 



5. Search for possible Lorentz invariance 
violation contributions. 



Two flavor v.. 



mass-induced oscillations 



are strongly favored over a wide range of alterna- 
tive explanations of the atmospheric v data |15) 
|16| . In a search for possible Lorentz invariance 
violation (LIV) contributions, we assumed mass- 
induced v oscillations as the leading mechanism 
and LIV as a sub-dominant effect I n this 

scenario one considers 2 flavor eigenstates, 2 mass 
eigenstates and 2 velocity eigenstates (character- 
ized by different maximum attainable velocities 
in the limit of infinite momentum) . In a first anal- 
ysis we considered the subsample of upthrough- 
going ju's for which the energy was estimated via 
MCS and we used the ratio R' 2 = N( ow /N' high , 
where low and high are for events with recon- 
structed energies E r v ec < 30 and > 130 GeV (av- 
erage energies of 13 and 146 GeV, respectively). 
In the analysis we fixed the neutrino mass os- 
cillation parameters at Am 2 = 0.0023 eV 2 and 
sin 2 2i? m = 1 and performed a minimization with 



respect to the LIV parameters Av and d v . The 
inclusion of LIV effects docs not improve the x 2 
of the fit and we obtain the 90% C.L. limits shown 
as the continuous line in Fig. [SJ The calculation 
was repeated for median neutrino energies of 17 
and 167 GeV, dashed line in Fig. [SJ 

We would like to acknowledge the cooperation 
of the members of the MACRO collaboration and 
discussions with several theoretical colleagues. 
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